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ABSTRACT 

In a quantum theory of gravity the gravitational Wilson loop, defined as a suitable quantum 
average of a parallel transport operator around a large near-planar loop, provides important in- 
formation about the large-scale curvature properties of the geometry. Here we shows that such 
properties can be systematically computed in the strong coupling limit of lattice regularized quan- 
tum gravity, by performing local averages over loop bivectors, and over lattice rotations, using an 
assumed near-uniform measure in group space. We then relate the resulting quantum averages to an 
expected semi-classical form valid for macroscopic observers, which leads to an identification of the 
gravitational correlation length appearing in the Wilson loop with an observed large-scale curva- 
ture. Our results suggest that strongly coupled gravity leads to a positively curved (De Sitter-like) 
quantum ground state, implying a positive effective cosmological constant at large distances. 
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1 Introduction 



An important question for any theory of quantum gravity is what gravitational observables should 
look like [1], i.e. which expectation values of operators (or ratios thereof) have meaning and physical 
interpretation in the context of a manifestly covariant formulation, specifically in a situation where 
metric fluctuations are not necessarily bounded. Such averages naturally include expectation values 
of the (integrated) scalar curvature and other related quantities (involving for example curvature- 
squared terms), as well as correlations of operators at fixed geodesic distance, sometimes referred 
to as bi-local operators. Another set of physical observables corresponds to the gravitational analog 
of the Wilson loop [2], providing information about the parallel transport of vectors, and therefore 
on the effective curvature, around large, near-planar loops, and the correlation between particle 
world-lines [3,4] (providing information about the static gravitational potential). In this paper we 
will concentrate on defining and exploring physical properties of the gravitational Wilson loop [5]. 

Before embarking on the gravitational case, let us recall more generally the well-known fact that 
many low energy physical properties in gauge theories cannot be computed reliably in weak coupling 
perturbation theory. Thus, for example, in non-Abelian gauge theories a confining potential for 
static sources placed in the fundamental representation is found at sufficiently strong coupling, by 
examining the behavior of the Wilson loop [2], defined for a large closed planar loop C as 



with = taA® and the t a 's the group generators of SU(N) in the fundamental representation. 
Specifically, in the pure gauge theory at strong coupling the leading contribution to the Wilson 
loop can be shown to follow an area law for sufficiently large loops 



where A(C) is the minimal area spanned by the planar loop C [6,7]. The quantity £ is the gauge 
field correlation length, defined for example from the exponential decay of the Euclidean correlation 
function of two infinitesimal loops separated by a distance |x|, 



Here the C e 's are two infinitesimal loops centered around x and respectively, suitably defined on 
the lattice as elementary square loops, and for which one has at sufficiently large separations 
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The inverse of the correlation length £ is known to correspond to the lowest mass excitation in the 
gauge theory, the scalar glueball. 

Not only will we adapt this definition to the gravitational case, but more specifically to the case 
of discrete gravity, in the context of the discretization scheme known as Regge calculus [8] . It will 
turn out that it is most easily achieved by using a slight variant of Regge calculus, in which the 
action coincides with the usual Regge action in the near-flat limit. In Section 2, we shall describe 
the lattice notion of parallel transport, and how areas are defined on the dual lattice. Then in 
Section 3, the gravitational Wilson loop will be defined, and it will be shown why, in the discrete 
case, we modify the action. In section 4, sample calculations will be performed and the behavior 
of large Wilson loops derived. Much of what is done will be in close parallel with the procedure in 
lattice gauge theories (as will be immediately obvious to those familiar with that area). We close 
with a discussion of the interpretation in semiclassical terms of our main result for a large Wilson 
loop. We will argue there that our results imply that for strong coupling (large bare Newton's 
constant G) the behavior of the Wilson loop is consistent with a positive vacuum curvature, and 
therefore with (Euclidean) De Sitter space. 



2 Rotations, parallel transport and Voronoi loops 



In lattice gravity, space-time is built up from flat simplices, with curvature restricted to the sub- 
spaces of codimension 2. In four dimensions, this means that the hinges, where the curvature lies, 
are triangles. The contribution to the action of a hinge is the product of its area, Ah, with the 
deficit angle, 5h , there. This is defined to be 2ir minus the sum of the dihedral angles at that hinge, 
in the simplices meeting there. We may also define a volume associated with each hinge, Vh (see 
later in this section). Then the lattice action for pure four-dimensional Euclidean gravity with a 
cosmological constant and the usual Einstein scalar curvature term is 

hatt = A J2V h - kJ2^hA h , (5) 
h h 

with k = l/(87rG), and Vh, Ah and Sh are all functions of the edge lengths, which are the basic 
variables in the theory, analogous to the metric in the continuum. This action only couples edges 
which belong either to the same simplex or to a set of neighboring simplices, and can therefore be 
considered as local, just like the continuum action. It leads to the regularized lattice functional 
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integral [9] 
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where, as customary, the lattice ultraviolet cutoff is set equal to one (i.e. all length scales are 
measured in units of the lattice cutoff). The lattice partition function Zi att should then be compared 
to the continuum Euclidean Feynman path integral 

Zcont = J [dg^v] exp j-A J dx Jg + j^-^ J dx^/gR^ . (7) 

In practice the lattice functional integral Zi a tt should be regarded as a regularized form of the 
continuum Euclidean Feynman path integral. The latter will involve a functional measure over 
metrics g^ u (x), usually of the form 

[dg,u] = II W)Y' 2 II d 9^) , (8) 



where a is a parameter, constrained by the requirement of a non-singular measure to a > — {d + 1). 
For a = \(d — 4){d + \) one has De Witt's measure, while for a = —(d + 1) one recovers the 
original Misner measure. Here we will mostly be interested in the physical four-dimensional case, 
for which d = 4 and therefore a = in the De Witt measure; the specific form of the functional 
measure is not expected to play an important role in the following, except that it will assumed to 
be diffeomorphism invariant [5]. 

Furthermore, unless stated otherwise, it will convenient to include the cosmological constant 
term in the measure as well, since this contribution is ultralocal and contains no derivatives of the 
metric, giving rise to an effective strong coupling measure dfj,(l 2 ), 

dfi{l 2 ) = [dl 2 ]e- x ^ h Vh . (9) 

This last expression represents a fairly non-trivial quantity, both in view of the relative complexity 
of the expression for the volume of a simplex, and because of the generalized triangle inequality 
constraints already implicit in the definition of [dl 2 ]. 

The main assumption used here regarding this effective strong coupling measure will be the 
existence of a stable ground state with a well-defined average lattice spacing, as implied by direct 
numerical evaluations of the lattice integrals in four dimensions, at least for sufficiently strong 
coupling [5,10]. In the following the lattice measure will therefore be assumed to be a suitable 
discretization of the continuum functional measure, and therefore of the form 

//*oo _ 
w 2 ] = / n ^ s )y n 4 @ ^ ■ ( io ) 

s ij 



with a again a real parameter, and a function of the squared edge lengths, ensuring the validity 
of the triangle inequalities. 

At strong coupling the measure and cosmological constant terms form the dominant part of the 
functional integral, since the Einstein part of the action is vanishingly small in this limit. Yet, and 
in contrast to strongly coupled lattice Yang-Mills theories, the functional integral is still non-trivial 
to compute analytically in this limit, mainly due to the triangle inequality constraints. Therefore, in 
order to be able to derive some analytical estimates for correlation functions in the strong coupling 
limit, one needs still to develop some set of approximation methods, which will discussed below. 
These methods and their results can later be tested by numerical means, for example by integrating 
directly over edges, through the explicit lattice measure over edges given above. 

One approach that appears natural in the gravity context follows along the lines of what is 
normally done in gauge theories, namely an integration over compact group variables, using the 
invariant measure over the gauge group [2]. It is of this method that we wish to take advantage here, 
as we believe that it is well suited for gravity as well. In order to apply such a technique to gravity 
one needs (i) to formulate the lattice theory in such a way that group variables are separated and 
therefore appear explicitly; (ii) integrate over the group variables using an invariant measure; and 
(iii) approximate the relevant correlation functions in such a way that the group integration can be 
performed exactly, using for example mean field methods for the parts that appear less tractable. 
In such a program one is aided, as will be shown below, by the fact that in the strong coupling limit 
one is expanding about a well defined ground state, and that the measure and the interactions are 
local, coupling only lattice variable (edges or rotations) which are a few lattice spacings apart, thus 
excluding the appearance of long-range (power-like) correlations. 

The downside of such methods is that one is no longer evaluating the functional integral for 
quantum gravity exactly, even in the strong coupling limit; the upside is that one obtains a clear 
analytical estimate, which later can be in principle systematically tested by numerical methods 
(which are exact). 

In the gravity case the analogs of the gauge variables of Yang-Mills theories are given by the 
connection, so it is natural therefore to look for a first order formulation of gravity. In a first order 
formalism one writes for the Einstein-Hilbert pure gravity Lagrangian density 



C = 



1 
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For any spacetime manifold with an affine connection one has for the Ricci tensor 



Rfiiy — 9 ° ' R\ 



(12) 
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where 

R\vo = 9 V T ^ — + r^r^ — r^r^ . (13) 

Variation of the pure gravitational action requires that 

^ J d 4 xS[^gg^R^} = • (14) 

The variation of the first and second terms inside the square parentheses are trivial, and the 
variation of R^ v can be simplified by virtue of the Palatini identity 

5R\ ua = sr^-sr^ . (15) 

After integrating by parts one can then show that the term involving the variation of the connection 
T implies 

dx9^ ~ 9ja^ x - 9f3(7 r° x = (16) 

(normally just written as g^ u; \ = 0), and which can be inverted to give the usual relationship 
between the connection T and the metric g in Riemannian geometry, namely 

^v = ^g Xa (d fl g ua + d u g II(7 -d a g^) . (17) 

Equating to zero the coefficients of 5g^ u gives instead the ten components of Einstein's field equa- 
tions. The (well-known) conclusion therefore is that in the quantum theory one can safely consider 
functionally integrating separately over the affine connection and the metric, treated as independent 
variables, with the correct relationship between metric and connection arising then as a consequence 
of the dynamics. For the rest of this section we will follow a similar spirit, separating out explicitly 
in the lattice action the degrees of freedom corresponding to local rotations (the analogs of the Ps 
in the continuum), which we will find to be most conveniently described by orthogonal matrices R, 
implying a choice of preferred coordinate systems within the simplices. 

The next step is a discussion of the properties of local rotation matrices in the context of the 
lattice theory, and how these relate to the lattice gravitational action. Since in Regge calculus 
the interior of each simplex s is assumed to be flat, one can assign to it a Lorentz frame S(s). 
Furthermore inside s one can define a ti-component vector <fi(s) = (<^>o • • • 4>d-i)- Under a Lorentz 
transformation of S(s), described by the d x d matrix A(s) satisfying the usual relation for Lorentz 
transformation matrices A T r\ A = rj, with r] the flat metric, the vector cf)(s) will rotate to 

4>'(s) = A(s)<f>(s) . (18) 

The base edge vectors ef = Iq^s) themselves are of course an example of such a vector. 
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Next consider two ti-simplices, individually labeled by s and s', sharing a common face f(s, s') 
of dimensionality d — 1 [10]. It will be convenient to label the d edges residing in the common face 
/ by indices i,j = 1 . . . d. Within the first simplex s one can then assign a Lorentz frame and 
similarly within the second s' one can assign the frame S(s'). The \d{d — 1) edge vectors on the 
common interface f(s, s') (corresponding physically to the same edges, viewed from two different 
coordinate systems) are expected to be related to each other by a Lorentz rotation R, 

%{s>) = R^(s',s)l^(s) . (19) 

Under individual Lorentz rotations in s and s' one has of course a corresponding change in R, 
namely R — > A(s') R(s', s) A(s). In the Euclidean d-dimensional case R is an orthogonal matrix, 
an element of the group SO(d). 

In the absence of torsion, one can use the matrix R(s', s) to describes the parallel transport of 
any vector <fi from simplex s to a neighboring simplex s', 

r(s') = R» V (8',S)F(8) (20) 

R therefore describes a lattice version of the connection [11,12]. Indeed in the continuum such a 
rotation would be described by the matrix 

R» v = [e Vdx Y v (21) 

with r*„ the affine connection. The coordinate increment dx is interpreted as joining the center 
of s to the center of s', thereby intersecting the face f(s,s'). Note that it is possible to choose 
coordinates so that R(s, s') is the unit matrix for one pair of simplices, but it will not then be unity 
for all other pairs. 

One can consider a sequence of rotations along an arbitrary path P(si, . . . , s„+i) going through 
simplices si . . . s n +i, whose combined rotation matrix is given by 

R(P) = R(s n+1 ,s n )---R(s 2 , Sl ) (22) 

and which describes the parallel transport of an arbitrary vector from the interior of simplex si to 
the interior of simplex s n+ i, 

^K+i) = RMP)</>"( Sl ) (23) 

If the initial and final simplices s n +i and s± coincide, one obtains a closed path C(s\, . . . , s n ), for 
which the associated expectation value can be considered as the gravitational analog of the Wilson 
loop. Its combined rotation is given by 

R(C) = R{s 1 ,s n )---R(s 2 ,s 1 ) (24) 
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Under Lorentz transformations within each simplex Sj along the path one has a pairwise cancellation 
of the A(sj) matrices except at the endpoints, giving in the closed loop case 

R(C) -► A( Sl )R(C)A T (si) (25) 

Clearly the deviation of the matrix R(C) from unity is a measure of curvature. Also, the trace 
trR(C) is independent of the choice of Lorentz frames. 




Figure 1. Elementary polygonal path around a hinge (triangle) in four dimensions. The hinge ABC, 
contained in the simplex ABCDE, is encircled by the polygonal path H connecting the surrounding vertices 
which reside in the dual lattice. One such vertex is contained within the simplex ABCDE. 

Of particular interest is the elementary loop associated with the smallest non-trivial, segmented 
parallel transport path one can build on the lattice. One such polygonal path in four dimensions 
is shown in Figure 1. In general consider a (d — 2)-dimensional simplex (hinge) h, which will be 
shared by a certain number m of (i-simplices, sequentially labeled by s% . . . s m , and whose common 
faces f(si, s%) ... /( ) f(s m , si) will also contain the hinge h. In four dimensions several 

four-simplices will contain, and therefore encircle, a given triangle (hinge). In three dimensions the 
path will encircle an edge, while in two dimensions it will encircle a site. Thus for each hinge h 
there is a unique elementary closed path C/j for which one again can define the ordered product 

R(C h ) = R(s 1 ,s rn )---R(s 2 ,s 1 ) (26) 

The hinge h, being geometrically an object of dimension (d — 2), is naturally represented by a 
tensor of rank (d — 2), referred to a coordinate system in h: an edge vector l 1 ^ in d = 3, and an area 
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bivector \{^h ~ tfj-h) in ^ = 4 etc. It will therefore be convenient to define a hinge bivector U in 
any dimension as 

U^ v {h) = M e^ vaiad _ 2 . . . l(£-2) ' (^) 
normalized, by the choice of the constant Af, in such a way that U 11V U ,J ' U = 2. In four dimensions 

U^(h) = ^L e ^ a pl?l$ (28) 

where l±(h) and I2Q1) two independent edge vectors associated with the hinge h, and Ah the area 
of the hinge. 

An important aspect related to the rotation of an arbitrary vector, when parallel transported 
around a hinge h, is the fact that, due to the hinge's intrinsic orientation, only components of the 
vector in the plane perpendicular to the hinge are affected. Since the direction of the hinge h is 
specified locally by the bivector U^ v of Eq. (28), one can write for the loop rotation matrix R 



R^(C) = (e ou ) r v (29) 

where C is now the small polygonal loop entangling the hinge h, and 5 the deficit angle at h. One 
particularly noteworthy aspect of this last result is the fact that the area of the loop C does not 
enter in the expression for the rotation matrix, only the deficit angle and the hinge direction. Note 
that in the above expression for the rotation matrix R(C) both the deficit angle 5(C) giving the 
magnitude of the rotation, as well as the bivector U (C) giving the direction of the rotation, are 
each rather complicated functions of the original edge lengths, with the latter also depending on a 
choice for the local coordinate system. 

At the same time, in the continuum a vector V carried around an infinitesimal loop of area Ac 
will change by 

AV^ = \R^ uXa A Xa V v (30) 

where A Xa is an area bivector in the plane of C, with squared magnitude A\ a A Xa = 2A C . Since 
the change in the vector V is given by SV a = (R — 1)°^ V 13 one is led to the identification 

\R a ^ v A^ = (R-l)^ . (31) 

Consequently the above change in V can equivalently be re-written in terms of the infinitesimal 
rotation matrix 

R^iC) = [e\ R ' A y (32) 
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(where the Riemann tensor appearing in the exponent on the r.h.s. should not be confused with 
the rotation matrix R on the l.h.s.). 

The area Ac is most suitably defined by introducing the notion of a dual lattice, i.e. a lattice 
constructed by assigning centers to the simplices, with the polygonal curve C connecting these 
centers sequentially, and then assigning an area to the interior of this curve. One possible way of 
assigning such centers is by introducing perpendicular bisectors to the faces of a simplex, and locat- 
ing the vertices of the dual lattice at their common intersection, a construction originally discussed 
in [13,14]. Another, and perhaps even simpler, possibility is to use a barycentric subdivision. Then 
the volume element, Vh, is defined by first joining the vertices of the polyhedron C, whose vertices 
lie in the dual lattice, with the vertices of the hinge h, and then computing its volume. We can then 
show that the polygonal area Ac is given by Ac{h) = d Vh/vjf 2 \ where vjf ^ is the volume of 
the hinge (a triangle area in four dimensions). 



3 Gravitational Wilson loop 



We have seen that with each neighboring pair of simplices s, s + 1 one can associate a Lorentz 
transformation R^ u (s, s + 1), which describes how a given vector transforms between the local 
coordinate systems in these two simplices, and that the above transformation is directly related to 
the continuum path-ordered (P) exponential of the integral of the local affine connection T* v (x) 
via 



m, = 



J path 
P g between simplices 



r x dx x 



(33) 



with the connection having support only on the common interface between the two simplices. Also, 
for a closed elementary path Ch encircling a hinge h and passing through each of the simplices that 
meet at that hinge one has for the total rotation matrix R = Y\ s R s , s +i associated with the given 
hinge 



[n 



R 



■s,s+l 



J(h)U(h) 



(34) 



as in Eq. (29). This matrix describes the parallel transport of a vector round the loop. 

More generally one might want to consider a near-planar, but non-infinitesimal, closed loop C, 
as shown in Figure 2. Along this closed loop the overall rotation matrix will still be given by 



RMC) = [n R 



s,s+l 



(35) 



sec 
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In analogy with the infinitesimal loop case, one would like to state that for the overall rotation 
matrix one has 



J{C)U{C)) 



(36) 



where £7 Mi ,(C) is now an area bivector perpendicular to the loop, which will work only if the loop 
is close to planar so that U^ u can be taken to be approximately constant along the path C . By 
a near-planar loop around the point P, we mean one that is constructed by drawing outgoing 
geodesies, on a plane through P. 




Figure 2. Gravitational analog of the Wilson loop. A vector is parallel-transported along the larger outer 
loop. The enclosed minimal surface is tiled with parallel transport polygons, here chosen to be triangles for 
illustrative purposes. For each link of the dual lattice, the elementary parallel transport matrices R(s, s') are 
represented by arrows. In spite of the fact that the (Lorentz) matrices R can fluctuate strongly in accordance 
with the local geometry, two contiguous, oppositely oriented arrows always give RRr 1 = 1. 

If that is true, then one can define an appropriate coordinate scalar by contracting the above 
rotation matrix R(C) with the some appropriate bivector, namely 

W{C) = u; a(3 (C)R at3 (C) (37) 

where the bivector, u a g(C), is intended as being representative of the overall geometric features of 
the loop. 

In the quantum theory one is of course interested in the average of the above loop operator 
W(C), as in Eq. (1). The previous construction is indeed quite analogous to the Wilson loop 
definition in ordinary lattice gauge theories [2], where it is defined via the trace of path ordered 
products of SU(N) color rotation matrices. In gravity though the Wilson loop does not give any 
information about the static potential [15,16]. It seems that the Wilson loop in gravity provides 
instead some insight into the large-scale curvature of the manifold, just as the infinitesimal loop 
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contribution entering the lattice action of Eqs. (5) and (39) provides, through its averages, insight 
into the very short distance, local curvature. 

Of course for any continuum manifold one can define locally the parallel transport of a vector 
around a near-planar loop C. Indeed parallel transporting a vector around a closed loop represents 
a suitable operational way of detecting curvature locally. If the curvature of the manifold is small, 
one can treat the larger loop the same way as the small one; then the expression of Eq. (36) for the 
rotation matrix R(C) associated with a near-planar loop can be re-written in terms of a surface 
integral of the large-scale Riemann tensor, projected along the surface area element bivector A a ^(C) 
associated with the loop, 



Thus a direct calculation of the Wilson loop provides a way of determining the effective curvature 
at large distance scales, even in the case where short distance fluctuations in the metric may be 
significant. Conversely, the rotation matrix appearing in the elementary Wilson loop of Eqs. (26) 
and (29) only provides information about the parallel transport of vectors around infinitesimal 
loops, with size comparable to the ultraviolet cutoff. 

Let us now look in detail at how to construct a Wilson loop in quantum gravity. Since this 
involves finding the expectation value of a product of rotation matrices round a loop, the natural 
procedure is to treat these rotation matrices as variables and to integrate over their product, 
weighted by the exponential of minus the Regge action. The expression for this action has been 
given in terms of functions of the edge lengths, but an alternative [11,12,17] is to find an expression 
for it in terms of the rotation matrices. For the dual loop around each hinge, the product of the 
rotation matrices gives the exponential of the deficit angle, 5, times the rotation generator, U, (see 
Eq. (29)) and we need to find a way of extracting the deficit angle from this product of matrices, 
at the same time as constructing a scalar function to be averaged. The obvious way of doing this 
is to contract the product of the i?-matrices with the rotation generator, U, and then take the 
trace. This is equivalent to the action used in [11,12,17] (see also [18,19] for a hypercubic lattice 
formulation), obtained by contracting the elementary rotation matrix R(C) of Eq. (29), with the 
hinge bivector of Eq. (27), 

U« 2 )=-^ E A h U aP {h)B^{h) (39) 

hinges h 

The above construction can be regarded as analogous to Wilson's lattice gauge theory, for which 
the action also involves traces of products of SU(N) color rotation matrices [2]. This contraction 




(38) 
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produces the sine of the deficit angle times the area of the triangular hinge and so for small deficit 
angles it is equivalent to the Regge action. However, in general, away from a situation of small 
curvatures, the two lattice action are not equivalent, as can be seen already in two dimensions. 

At this stage, we choose to differ from the choice made in [17], for the simple reason that when 
we come to evaluate Wilson loops, the final result often involves the trace of the bivector U, which 
is zero. Therefore, instead of contracting with U, we use a linear combination of it and the unit 
matrix. In particular, we take the contribution to the action, of a hinge, labelled h, to be 

h = \ A h Tr[(U h + e h) (R h - R^ 1 )], (40) 

where e is an arbitrary multiple of the unit matrix in four dimensions. We have subtracted the 
inverse of the rotation matrix for the hinge for reasons that will become apparent when we evaluate 
Wilson loops, and it also plays an important role in the action. At the end, we shall be interested 
in the limit of small but non-zero e. 

The classical action may be evaluated as follows. Since R equals the exponential of 8 times U, 
it may be expanded in a power series in 5, which is then contracted with the (U + e/4) and the 
trace taken. We use 

Tr(U 2n+1 ) = ; Tr(U 2n ) = 2 (-l) n , (41) 

to show that 

I h = - k A h sin(^), (42) 

independently of the value of the parameter e. (The unit matrix times the even terms in the 
power series expansion produces a cos(5) term, but this cancels between the R and the R -1 
contributions.) Thus e is in fact an arbitrary parameter, which can be conveniently taken to be 
non-zero, as we shall see. 

There is now a slight amount of freedom in how we define the Wilson loop, for a path C in the 
dual lattice of a simplicial space. The main choices seem to be 

(i) W(C) = < Tr{Ri R 2 ... R n ) >; (43) 
(ii) W(C) = < Tr[(U c + e J 4 ) Ry R 2 Rn] > ■ (44) 
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Here the Ri are the rotation matrices along the path; in (ii), there is a factor of (Uc + CI4), 
containing some "average" direction bivector, Uc, for the loop, which, after all, is assumed to be 
almost planar. The position of the Uc term in the product of R^s is not arbitrary; to give a unique 
answer, it needs to be placed before an R which begins one of the plaquette contributions to the 
action. 

We would like to take as independent fluctuating variables the rotation matrices Ri and the loop 
bivectors (7j, in a first order formalism similar in spirit to that used in [17]. This last statement 
clearly requires some clarification, as both the rotation matrices and the loop bivectors depend 
on the choice of the original edge lengths, as well as on the orientation of the local coordinate 
system, and cannot therefore in general be considered as independent variables (as should have 
already been clear from the detailed discussion of the properties of rotation matrices given in the 
previous section). On the lattice strong edge length fluctuations get reflected in large fluctuations 
in the local geometry, which in turn imply large correlated fluctuations in both the deficit angles 
and in the orientations of the elementary loop. It would therefore seem at first that one would 
have to integrate over both sets of coupled variables simultaneously, with some non-trivial measure 
derived from the original lattice measure over edge lengths, which in turn would make the problem 
of computing the Wilson loop close to intractable, even in the strong coupling limit. In particular 
one has to take notice of the fact that the lattice deficit angles and the loop bivectors are related to 
the metric and connection, as they appear in a first order formulation, in a rather non-trivial way. 

But there are two important aspect that come into play when evaluating the expectation value 
of the gravitational Wilson loop for strongly coupled gravity, the first one being that the overall 
geometric features of the large near-planar loop provide a natural orientation, specified for example 
by a global loop bivector Uc- As will become clear from explicit calculations given below, in the 
strong coupling limit the tiling of the large Wilson loop surface by elementary parallel transport 
loops, which in general have random orientations, requires that their normals be preferentially 
oriented perpendicular to the plane of the loop, since otherwise a non-minimal surface must result, 
which leads to a necessarily higher order contribution in the strong coupling limit. 

In the case of a hinge surrounded by the large loop with bivector Uc, one is therefore allowed 
to write for the bivector operator Uh associated with that hinge, labelled by h, 

U h = U c + SU h (45) 

where SUh is the quantum fluctuation associated with hinge bivector at h. But assuming the 
fluctuation in SUh to be zero is an unnecessarily strong requirement, and in the following it will 
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be sufficient to take < 5Uh >= and < (SUh) 2 >^ 0, which can be regarded as a mean-field type 
treatment for the loop bivectors. It will be important therefore in the following to keep in mind 
this distinction between the fluctutating hinge bivector Uh, and its quantum average. 

The second important aspect of the calculation is that at strong coupling the edge lengths, 
and therefore the local geometry, fluctuate in a way that is uncorrelated over distances greater 
than a few lattice spacing. Thus, mainly due to the ultralocal nature of the gravitational lattice 
measure at strong coupling, the fluctuations in the U's can be taken as essentially uncorrelated as 
well, again over distances greater than a few lattice spacings, which further simplifies the problem 
considerably. 

One would expect that for a geometry fluctuating strongly at short distances (corresponding 
therefore to the small k limit) the infinitesimal parallel transport matrices R(s, s') should be dis- 
tributed close to randomly, with a measure close to the uniform Haar measure, and with little 
correlation between neighboring hinges. In such instance one would have for the local quantum 
averages of the infinitesimal lattice parallel transports < R >= 0, but < R R^ 1 0, which would 
require, for a non-vanishing lowest order contribution to the Wilson loop, that the loop at least be 
tiled by elementary loops with action contributions from Eqs. (5) or (39), thus forming a minimal 
surface spanning the loop. Then, in close analogy to the Yang-Mills case of Eq. (2) (as a general 
reference, see for example [20]), the leading contribution to the gravitational Wilson loop would be 
expected to follow an area law, 

<W{C)>~ const. k Mc ^ ~ exp(-A(C)/£ 2 ) (46) 

where A(C)) is the minimal physical area spanned by the near-planar loop C, and £ the gravitational 



correlation length, equal to £ = In fc| for small k. For a close-to-circular loop of perimeter P 
one would use A(C) « P 2 /4tt. 

We choose now to focus on the Euclidean case in four dimensions, where the rotation matrices 
will be elements of SO (A). In evaluating the averages over the rotation matrices in the expectation 
values in the Wilson loops, the integrations we have to perform will be of the form 

J (f[dfi H (Ri)j 2V[...(^ + eJ 4 )...i2 fc ...] exp I -| £ A h Tr[(U h + eh)(R h - R^ 1 )] j / J\f, 

\i=l / y hinges h J 

(47) 

where the normalization factor is given by 
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AT = J (fldfiHiRiU expj-^ £ A h Tr[(U h + eI 4 )(K h --R^ 1 )} 

\i=l / \ hinges h 



(48) 



hinges h 

This factor will be omitted from subsequent expressions, for notational simplicity. 

For smooth enough geometries, with small curvatures, the rotation matrices can be chosen to 
be close to the identity. Small fluctuations in the geometry will then imply small deviations in 
the i?'s from the identity matrix. However, for strong coupling (k — ► 0) the usual lattice measure 
/ d/j,(l 2 ) [21] does not significantly restrict fluctuations in the lattice metric field. As a result we 
will assume that these fields can be regarded, at least in this regime, as basically unconstrained 
random variables, only subject to the relatively mild constraints implicit in the measure dfi(l 2 ). 
Thus as k — > 0, the geometry is generally far from smooth, since there is no coupling term to 
enforce long range order (the coefficient of the lattice Einstein term goes to zero) , and one has as a 
consequence large local fluctuations in the geometry. The matrices R will therefore fluctuate with 
the local geometry, and average out to zero, or a value close to zero, in the sense that, for example, 
the 50 (4) rotation 



Re 



/cos 9 -sinO 0\ 

sin 9 cos 9 

10 

V 1/ 



(49) 



averages out to zero when integrated over 9. In general an element of SO(n) is described by 
n(n — l)/2 independent parameters, which in the case at hand can be conveniently chosen as the 
six 50(4) Euler angles. The uniform (Haar) measure over the group is then 

j /*27T /*7T f'7T f'7T f'7T f-7l 

dn H {R) = □ / Mi d9 2 d9 3 d9 4 sin 9 4 d9 5 sin 9 5 d9 6 sin 2 9 6 (50) 

327r y Jo Jo Jo Jo Jo Jo 

This is just a special case of the general n result, which reads 

(n \ n—1 i 

n r(i/2)/2 n ^(«+ 1 )/ 2 n n sw?- 1 ^ ^ (51) 
i=i / i=i j=i 

with < 9\ < 2vr, < 9 j k < ir [22]. 

These averaging properties of rotations are quite similar to what happens in SU(N) Yang- 
Mills theories, or even more simply in (compact) QED, where the analogs of the 50(d) rotation 
matrices R are phase factors U^x) = e iaA ^ x \ There one has the trivial group averaging property 
/ Ufj,(x) = and / U^x) UUx) = 1. In addition, for two contiguous closed paths C\ and 
Ci sharing a common side one has 

e ^ Cl A - dl e ^C 2 A - dl = e ^ c A-dl = e ,/ s B.ndA ^ (52) 
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with C the slightly larger path encircling the two loops. For a closed surface tiled with many 
contiguous infinitesimal closed loops the last expression evaluates to 1, due to the divergence 
theorem. In the lattice gravity case the discrete analog of this last result represents the (exact) 
lattice analog of the contracted Bianchi identities [23] . 

In practice, luckily, we do not have to explicitly integrate over SO (4) angles, but rather use the 
following properties of the Haar measure, normalized to one, on the group: 



/ 



d H R = 1 ; (53) 



J d H R Tr(A R) TriR- 1 B) = - Tr(A B) , (54) 
for arbitrary 4x4 matrices A and B, which also implies 

J d H R Rij Rkl = ^ Su S jk . (55) 

As stated previously, we will regard the individual hinge bivectors Uh as aligned on average 
with the Wilson loop bivector Uc- Alternatively, one can perform the following simple exercise, in 
which one assumes each hinge bivector Uh is aligned in general (arbitrary) directions, and performs 
an integration over those directions. Here we will given an example of such a calculation. 

Of course the normals to a 2-d loop in a 4-d space form a plane, rather than a single direction, 
so in general what one needs to do is an integration over the directions in that plane. However in 
this example, we are not restricting the plane of the loop, so the integration is a four-dimensional 
one over all possible directions. It is performed as follows: in the definition of U a p in Eq. (28), we 
put 

lf = \h\a a , l$ = \l 2 \bP, (56) 

with a 7 a 7 = 6 7 6 7 = 1. Since Ah = ^iH^I sin<^>, where <j) is the angle between the vectors l\ and 
I2, the dependence on the magnitudes of those two vectors cancels from the expression for U and 
we obtain 

U a p = 7 — , (57) 
sm0 

with cos0 = a 7 6 7 . The integration of U a p then becomes 
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i n /' *, n /' ^w^w-D^-D, (58) 



i=l 

where the coefficient in front of the integral ensures that the measure is normalized to 1. 



4 Evaluation of Wilson loops and behavior for large loops 

In this section, we first evaluate < W > for some simple loops and then discuss the general behavior 
for arbitrary loops, ending with a consideration of the asymptotic behavior for large loops. We 
shall work out < W > for the two possible definitions listed in the previous section. 

4.1 Loop round a single hinge 

Consider a single hinge of area A, at which four 4-simplices meet (see Figure 3.). The loop C will 
consist of four segments between the Voronoi centers of the simplices. Let the rotation matrices on 
these segments be R\, R 2 , R3, i?4, and the rotation generator for the hinge U. Then we have either 



or 



(i) W(C) = < Tr(R\ R 2 R 3 R4) > 



(it) W(C) =< Tr[(U + e I 4 ) Ri R 2 R3 Ri] > 



(59) 
(60) 




Figure 3. A parallel transport loop with four oriented links on the boundary. The parallel transport matrices 
R along the links, represented here by arrows, appear in pairs and are sequentially integrated over using the 
uniform measure. 
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Since the only non-vanishing contribution to the integration over the i?'s will come from the 
product of an Ri with the corresponding R^ 1 , then the lowest order contribution in k will come 
from the term in the expansion of the exponential of minus the action which is linear in R _1 . Thus 
in case (i) we obtain 

k f 

-A d H R x d H R 2 d H R 3 d H R A Tr(Ri R 2 R 3 R A ) Tr[(U + e h) R^ 1 R^ 1 R 2 l R[\ (61) 

and a similar expression for (ii), with the extra factor of (U + e I4) inserted. Integration over the 
i?'s results in 

(») \\ATr(U + eh) = \Ae (62) 

or 

(ii) \\ATr{U 2 + 6 2 I 4 ) = -\A{1 - 2e 2 ), (63) 
and integration over the C/'s (a sum over all possible orientations of the loop) is trivial. 

4.2 Loop round two hinges 

Suppose now that the loop goes around two adjacent hinges, with rotation bivectors Ui,U2- The 
R matrices going clockwise round the loop in Figure 4., starting at the top left hand corner, are 
labelled 1, 2, 5, 6, 7, 4 and the one between the two loops is R3. 




Figure 4. A parallel transport loop with six oriented links on the boundary. 
Then the possible values of the Wilson loop are 



(i) W{C) = < Tr(R 1 R 2 R 5 R 6 R 7 R4) >, (64) 



19 



or 



(U) W{C) = < Tr[{U c + e J 4 ) R-y R 2 R 5 Re R7 Ri] > • (65) 

This time, the lowest order non-zero contribution will come from the term in the expansion of 
the exponential which is quadratic and involves the product of the R _1 's from the two hinges, so 
for (i) the integration over the i?'s gives 

k 2 f 

— A 1 A 2 j d H R 1 dHR2dHRzdHRid H R^d H R G d H R 1 Tr{R 1 R 2 R^R Q R 7 R i ) 
x Tr[([/i + e h) i?7 1 R^ 1 R^ 1 R^ 1 ] Tr[(U 2 + e h) R 3 Rj 1 Rq 1 R^ 1 ] , (66) 

and a similar expression for (ii). Evaluation of the integrals leads to 

(i) ^^A 1 A 2 Tr(lq + e h) Tr(U 2 + e J 4 ) = ^ A 1 A 2 e 2 , (67) 

or 



(**) ^AiA 2 Tr[(U c + e / 4 ) (I7i + e / 4 )] Tr(f7 2 + e J 4 ) = ^ Ai A 2 e [Tr(C/ c C/i) + 4 e 2 ] . 

(68) 

The integration over the ?7's, and therefore over the loop's orientation, in Eq. (67) is trivial. 
In Eq. (68), we either set U\ = Uc + SU±, with < 5U\ >= after which the sum over the loop's 
orientation also becomes trivial since Tr(UcU\) = —2, or we can use the following technique to 
evaluate the integral of Tr(UcU\): we assume (as usual) that Uc corresponds to a planar or almost- 
planar loop and we choose Cartesian coordinates (in the tangent space to it) such that the loop is 
in the (1,2) plane. We take the generator of Uc to be a right-angled triangle in the (3,4) plane, 
with l C \ = (0, 0, 1, 0), l C2 = (0, 0, 0, 1) and A c = 1/2. It can then be shown that 

Tr{U c U{) = ~ [(h)3(l2)4 ~ (h)4(h)3\ , (69) 

where l\ and l 2 are edge- vectors for the triangular hinge for U\ and A is the triangle's area. The 
integral to be evaluated here is therefore 

A n r n r 1)^-1). (70) 
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It can be shown that this integral is zero (an indication that this might be so is the antisymmetry 
in a and b in the integrand), and so Eq. (68) becomes identical to Eq. (67), apart from an extra 
power of e. 

4.3 Loop with one internal hinge 

We now need to consider the situation where the Wilson loop goes around a number of hinges and 
there is at least one internal hinge, i.e. a hinge where the elementary loop surrounding it is not part 
of the Wilson loop. For simplicity, we shall consider the case of one such loop. For the labelling of 
the rotation matrices and the hinges, the reader can annotate Figure 5. in a way consistent with 
the expressions below. 




Figure 5. A larger parallel transport loop with twelve oriented links on the boundary. As before, the parallel 
transport matrices along the links appear in pairs and are sequentially integrated over using the uniform 
measure. The new ingredient in this configuration is an elementary loop at the center not touching the 
boundary. 

In this case, the lowest order contribution comes from a ninth-order term in the expansion of 
the exponential of the action. We obtain the following results in the two cases: 

W J 4^ (A A ) (ft Tr(U t + el 4 )) = jl (j[ A^j <?, (71) 
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or 



* 9 1 

(") 49 4^ 



' 9 \ / 9 N 

IJ ^ rr[(l7 c + e/4) (l/i + eh)} II Tr (^ + e/ 4) 
vi=l / \i=2 y 



A "' '[[4, e 8 [rr(C/ c C/ 1 )+4e 2 ]. (72) 
i=i / 
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Integration over the U's is trivial for Eq. (71). In Eq. (68), we either set U\ = Uc + SUi, 
with < 6U± >= after which the sum over the loop's orientation also becomes trivial, or we can 
integrate over the relative orientation of U\ relative to Uc using the integral formula given at the 
end of the previous section, after which Eq. (72) gives the same as Eq. (71) but for an extra power 
of e. The above result also shows that it is better to take e > 0, otherwise the answer vanishes to 
this order. But this is not a problem, as the correct lattice action is recovered irrespective of the 
value of e, as shown earlier in Eq. (42). 



4.4 Large loop 

The value of a Wilson loop, in the case when the loop is very large and surrounds n hinges, can be 
seen to be roughly of the general form 

ii(n^)^b + ^r, (73) 

where a + (3 = n. If A is of the order of the geometric or arithmetic mean of the individual loops, 
this can be approximated by 

(^4) n ea[p+q ' 2]P - (74) 

The above result shows again that one should take e > 0, otherwise the answer vanishes to this 
order. As mentioned previously this is quite legitimate, as the correct lattice action is recovered 
irrespective of the value of e, as in Eq. (42). Then using n = Ac /A, we may write the area- 
dependent first factor as 

exp[ (Ac /A) log(fc A/16) ] = exp(- A c /£ 2 ) (75) 
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where £ = y[A/\ log(k A/16)\]. Recall that this is in the case of strong coupling, when k — > 0. The 
above is the main result of this paper. The rapid decay of the quantum gravitational Wilson loop as 
a function of the area is seen here simply as a general and direct consequence of the assumed disorder 
in the uncorrelated fluctuations of the parallel transport matrices R(s, s') at strong coupling. 

We note here as well that the correlation length £ is defined independently of the expectation 
value of the Wilson loop. Indeed a key quantity in gauge theories as well as gravity is the correlation 
between different plaquettes, which in simplicial gravity is given by 



< (6A) h {6A) h , > 



d/i(l 2 )(8A) h (5A) h , e k ^ A " 



dfi(l 2 )e k ^ 5hAh 



(76) 



In order to achieve a non- vanishing correlation one needs, at least to lowest order, to connect the 
two hinges by a narrow tube [24], so that 



< {5A) h {5A) h , > c ~ (k nt ) 1 ~ e - rf (WC 



(77) 



where tit I represents the minimal number of dual lattice polygons needed to form a closed surface 
connecting the hinges h and h' (as an example, for a narrow tube made out of cubes connecting 
two squares one has rat =4). In the above expression d(h, h') represents the actual physical distance 
between the two hinges, and the correlation length is given in this limit (A; — > 0) by £ ~ lo/nt \ log k\. 
where Iq is the average lattice spacing. Here we have used the usual definition of the correlation 
length £, namely that a generic correlation function is expected to decay as exp(— distance/^) for 
large separations. Figure 6. provides an illustration of the situation. 




Figure 6. Correlations between action contributions on hinge h and hinge h! arise to lowest order in the 
strong coupling expansions from diagrams describing a narrow tube connecting the two hinges. Here vertices 
represent points in the dual lattice, with the tube-like closed surface tiled with parallel transport polygons. 
For each link of the dual lattice, the 50(4) parallel transport matrices R are represented by an arrow. 
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The strong coupling area law behavior predicted for a large Wilson loop in Eq. (75) should 
be compared with the results for this in numerical simulations of lattice gravity. For small deficit 
angles (small curvature), the action used in this paper [involving Eq. (40)] is sufficiently close 
to the usual Regge action of Eq. (5) that the standard simulations can be used for comparison. 
Universality arguments would suggest a similar behavior for the gravitational Wilson loop for a 
wide class of lattice actions, constructed so as to reproduce the Einstein-Hilbert continuum action 
in the continuum limit. 



5 Conclusions 



In this paper, we have provided two possible constructions of Wilson loops in the gravitational 
case, and have given explicit calculations for small loops, as well as deriving the asymptotic be- 
havior for large loops. The final step, which we will give here, will be an attempt at providing an 
interpretation of this last and main result in semiclassical terms. As discussed in the introduction, 
the rotation matrix appearing in the gravitational Wilson loop can be related classically to a well- 
defined physical process: a vector is parallel transported around a large loop, and at the end it is 
compared to its original orientation. The vector's rotation is then directly related to some sort of 
average curvature enclosed by the loop. The total rotation matrix R(C) is given in general by a 
path-ordered (V) exponential of the integral of the affine connection T* v via 



R a JC) = \Vexp\(f r x .dx x \} a . (78) 



Upath C 



In such a semiclassical description of the parallel transport process of a vector around a very large 
loop, one can re-express the connection in terms of a suitable coarse-grained, or semi-classical, 
Riemann tensor, using Stokes' theorem 



I JS(C) 



R%{C) ~ [exp 1 / _ R-.^Ag T , (79) 



where here A^S is the usual area bivector associated with the loop in question, 

A% = \ I dx" x v '. (80) 



The use of semi-classical arguments in relating the above rotation matrix R(C) to the surface 
integral of the Riemann tensor assumes (as usual in the classical context) that the curvature is 
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slowly varying on the scale of the very large loop. Since the rotation is small for weak curvatures, 
one can write 



i?yc) ~ [1 + 1 f R; lw Ag + ...] a fj . (81) 



IS(C) 

At this stage one is ready to compare the above expression to the quantum result of Eq. (75), and 
in particular one should relate the coefficients of the area terms, which leads to the identification 
of the magnitude of the large scale semiclassical curvature with the genuinely quantum quantity 
l/£ 2 . Since one expression [Eq. (81)] is a matrix and the other [Eq. (75)] is a scalar, we shall take 
the trace after first contracting the rotation matrix with (Uq + e/4), as in our second definition of 
the Wilson loop, giving 



W(C) ~ Tr ([U c + eh) exp j \ jf 



*:„*S • <«) 



As is standard in simplicial gravity, we use 



Ra j3fiu — RU a pU^ v , (83) 

where R is some average curvature over the loop, and the f7's here will be taken to coincide with 
Uc- The trace of the product of (Uc + with this expression gives 

Tr(RU c A c ) = - 2RA C , (84) 

where we have used U^ v A c v = 2 Ac (we choose the directions of the bivectors such that the latter 
is true for all loops). This is to be compared with the linear term from the other exponential 
expression, — Ac/i 2 . Thus the average curvature is computed to be of the order R ~ l/£ 2 , at 
least in the small k limit [25] . An equivalent way of phrasing the last result is that l/£ 2 should be 
identified, up to a constant of proportionality, with the scaled cosmological constant A, which can 
be regarded as a measure of the intrinsic curvature of the vacuum. 

One important aspect of lattice gravity, and of the estimate for the large scale behavior of the 
Wilson loop given here, is that it can be tested by numerical methods. That is, the large scale 
behavior of the Wilson loop can in principle be computed directly by evaluating the lattice path 
integral using numerical methods, thus bypassing entirely the need for a separate treatment of the 
quantum fluctuations in the rotation matrices R and in the elementary loop bivectors U, as was 
done here in order to obtain an analytical result in the strong coupling limit. Furthermore in a 
numerical treatment of the Wilson loop one is no longer restricted necessarily to the strong coupling 
limit. 
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We see that a direct calculation of the Wilson loop for gravity can provide an insight into 
whether the manifold is De Sitter or anti-De Sitter at large distances 3 . Note that the definition of 
the gravitational Wilson loop is based on a surface with a given boundary C, in the simplest case 
the minimal surface spanning the loop. It is possible though to consider other surfaces built out of 
elementary parallel transport loops. These will be considered elsewhere. 
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3 In contrast, numerical studies of full lattice quantum gravity show that the elementary Wilson loop, describing 
the parallel transport around a single hinge, provides evidence for negative curvature at distance scales comparable 
to the ultraviolet cutoff [5] . 
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